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The structure of the Meissner effect in a current-carrying cylindrical wire with arbitrary disorder is studied 
following a numerical procedure that is exact within the quasiclassical approximation. A distribution of 
current is found that is non-monotonous as a function of the radial coordinate. For high currents, a robust 
gapless superconducting state develops at the surface of both clean and dirty leads. Our calculation provides 
a quantitative theory of the critical current in realistic wires. 
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The Meissner effect is one of the most fundamental 
properties of the superconducting state. It originates 
from the existence of a disipationless current density gen- 
erated by the superfluid velocity distribution which char- 
acterizes the state of broken gauge symmetry. The cur- 
rent density j depends on the velocity Vs through a com- 
plicated functional j[vs] which in general is non-linear 
and non-local The Pippard approximation assumes 
a linear relation j(r) = J K(r,r') • Vs{r')dr', which is 
valid for low enough superfluid velocities ||^ and which 
in the local limit reduces to London's equation. A local 
but yet non-linear approximation to the full functional 
j[vs] is implemented, for instance, within the context of 
a Ginzburg-Landau description . 

In this Letter, we present a fully non-linear and non- 
local numerical study of the field and current distribu- 
tions in a current-carrying cylindrical wire. Our calcula- 
tion of the hdl functional j[vs] is exact within the frame- 
work of the quasiclassical approximation. We encounter 
a rich physical structure determined by the interplay be- 
tween non-linearity, non-locality, and the global stability 
of the current configuration. In particular we find that, 
as a function of the distance to the central axis, the cur- 
rent density is non-monotonous for currents close to the 
critical value, displaying a maximum near the surface. 
This configuration can be viewed as precursor of the in- 
termediate state j^. We also find that, for high total 
currents, the superfluid velocity near the surface acquires 
values so large that they cannot be realized in a quasi- 
one-dimensional wire. In a three-dimensional wire, these 
large values of the superfluid velocity are possible because 
they are supported by the global stability of the current 
distribution. This causes a strong distortion of the lo- 
cal quasiparticle spectrum, which then develops a robust 
gapless form. The calculation presented here provides a 
quantitative theory of the critical current and represents 
an improvement over the phenomenological Silsbee's cri- 
terion [||. 



We investigate the structure of currents and fields in a 
cylindrical wire made of a type I s-wave superconductor 
with an arbitrary degree of disorder due to non-magnetic 
impurities. We have explored a broad range of temper- 
atures and wire radii, and have found that the most in- 
teresting physics appears at low temperatures and for 
large wire radii (specifically, for radius R > S,o-, Xq, where 
^0 and Aq are the zero temperature coherence and pen- 
etration lengths). Thus we have focussed on wires that 
are wide enough both to let the different physical mag- 
nitudes vary across its section, and to prevent thermal 
or quantum phase slips from taking place. We envis- 
age a steady-state scenario where, in the absence of ex- 
ternally applied fields, a disipationless current distribu- 
tion j[vs] (r) flows through the wire generating a magnetic 
field B(r) = {mc/e)\7 x v^, where the superfluid veloc- 
ity is Vs = {h/2m)[V(f> — {2e/ch)A]. Ampere's law then 
reads 
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At this point wc make use of the BCS theory of su- 
perconductivity in its quasi-classical formulation 
Together with Eq. (|l|), the following set of equations 
must be solved self-consistently: 
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Here, g(r,p, i?) is the quasiclassical fermion propagator 
in Keldysh space, whose components are g^, g^, and 
E = Sbcs + Simp is its self-energy, which has 
contributions Q from the pairing interaction (Sbcs) and 



1 



from impurities (Ejmp), p is a unit vector in the direction 
of ki?, A is the modulus of the superconducting order 
parameter, U is the coupHng constant, and A^o is the 
normal density of states, are Pauli matrices and fa 
is a block-diagonal matrix with block entries like in 73. 
Disorder is included through Simp and is characterized 
by the dimensionless parameter T = H/t Aq, t being the 
elastic scattering time and Aq the zero temperature, zero 
current superconducting gap. Equations (|l]), (^) and (|^), 
together with the continuity equation, can be derived as 
the time-independent equations satisfied by the extrema 
of the gauge-invariant action of an s-wave superconductor 
whose fermion propagator obeys Eq. (^). 

We exploit the cylindrical symmetry of the problem 
and consider configurations in which j points parallel to 
the wire axis. The continuity equation can then be shown 
to imply that is also directed along the wire, so that 
the magnetic field B has only angular component. Like 
the order parameter A, these three quantities depend 
only on the radial coordinate p. 

The self-consistency equations are solved according to 
the following procedure: We assume initial profiles A(p) 
and Vs{p) for the order parameter and the superfluid ve- 
locity. The quasiclassical fermion propagator is now ob- 
tained for all values of p,p,E by solving the equation 
of motion (|^) with hard-wall boundary conditions. This 
requires a self-consistent calculation of Simp(-B) for all 
values of E. The resulting is introduced into Eqs. 

and to determine A(p) and j{p)- The superfluid 
velocity Vg (p) is then obtained by solving the differential 
equation (|l|). The described calculation casts A(p) and 
Vs{p) as output distributions which are reintroduced as 
input for the next iterative step. The whole procedure is 
repeated until self-consistency in A and Vs is achieved. 
The value of Vs{p = 0) is given as an initial input param- 
eter and is kept fixed throughout the succesive iterative 
steps, effectively acting as a label for the resulting current 
configuration. 

Except where stated otherwise, we present results for 
a superconducting wire of radius i? = 5 Aq at a tem- 
perature T ~ 0.2 Tc- The Ginzburg-Landau parameter is 
K = Ao/^o = 0.4 and the Debye energy is Ejj = 10 Aq. In 
Fig. 1, we plot the total current / as a function of the su- 
perfluid velocity at three different locations in the wire 
(p = 0, i?/2, and R). Current densities and velocities are 
rather small in the core of the wire, even when / equals 
its critical value Ic- In contrast to this, Vs can be quite 
large near the surface. Fig. 1 indicates that it can amply 
exceed the stability threshold of a quasi-onc-dimensional 
wire (i? ^ Ao, Co) made of the same material. These large 
values of Vg are possible because it is the global stability 
of the current configuration what matters. Thus Vs can 
be very high in a small region of space, provided that it is 
sufficiently small in other regions in such a way that the 
configuration is globally stable. In Fig. 1 we plot I{vs) 
only for the clean (F = 0) case, but qualitatively similar 
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FIG. 1. Total current 7, in units of Iq = STreA'oWFAoAo, 
plotted as a function of the superfluid velocity Vs{p), in units 
of Vd = Ao/pf, at three diff'erent values of p, for a clean 
(r = 0) wire of radius _R = 5Ao at a temperature T = 0.2 Tc. 
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FIG. 2. (a) Order parameter A, (b) superfluid velocity «s, 
(c) magnetic field B, and (d) current density j, as a func- 
tion of the radial coordinate p, for the same wire as in Fig. 
1. The different current configurations are characterized by 
Ws(0). The units for current density and magnetic field are 
jo = 2eA'^o«FAo/3 and Bo = GneNovpAoXo/c. The bulk crit- 
ical field He ~ 0.54i3o is indicated in (c). Inset in (d); j{p) 
at 7 = 7c for several disorder strengths. 
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results are obtained for dirty wires. 

In Fig. 2 we present profiles of the order parameter, 
superfluid velocity, magnetic field, and current density, 
for several current configurations. For low currents, the 
order parameter A{p) stays essentially flat, since Vs is 
small everywhere. As the current increases, A develops 
a depression in the vicinity of the surface, brought about 
by the large values which Vg acquires in that region. Fig. 
2(b) shows a monotonous increase of Vs as a function of p. 
This is consistent with the constant sign of the magnetic 
field B shown in Fig. 2(c), since the two quantities are 
related by dvg/dp = {e/mc)B. In turn, B{p) also shows 
a monotonous increase which however tends to level off 
near the surface. Fig. 2(c) shows that, at the surface, 
B can exceed the thermodynamic critical field He by as 
much as a factor of 2, in clear violation of Silsbee's cri- 
terion, which states that the critical current is achieved 
when B{p = R) equals He. 

The large values which, for high current configurations, 
the superfluid velocity is forced to adopt near the surface 
are responsible for some of the most interesting physical 
features of the Meissner problem. The current density 
displays a peculiar behavior that can be approximately 
understood by noting that, in a quasi-one-dimensional 
wire, j is not a monotonous function of Vs but instead 
displays a maximum 10|. Thus, when Vs becomes 
very large, j no longer increases with Vg but actually 
decreases. Translated to the three-dimensional wire, the 
consequence is that, as a function of p, j goes through 
a maximum near the surface. This can be clearly ap- 
preciated in Fig. 2(d), whose inset shows in addition 
that a similar (but yet less pronounced) high-current be- 
havior can be displayed in the presence of appreciable 
disorder. We remark that this phenomenon is realized 
for local values of Vg that would be unstable (or directly 
not realizable) in a quasi-one-dimensional wire. 

The explanation given above suggests that the peak in 
the current density profile can be understood in terms 
of a strictly local theory which invokes the j{vs) relation 
of the quasi-onc-dimcnsional wire. This is the case, for 
instance, in a Ginzburg-Landau description, where the 
local relation j{vs) oc Vs(l — vf) holds. We can check 
the adequacy of the local hypothesis, by considering the 
ansatz j oc |ApUs, with A and Vs taken from the exact 
calculation. The result (not shown) is that, while the 
qualitatively correct behavior is reproduced in the dirty 
limit, the peak structure disappears in the clean case. 
The conclusion is that an adequate description of the 
peak structure in the current profile of a wire requires in 
general the combined inclusion of non-locality and non- 
linearity in an essential way. 

From the experimental viewpoint, the most relevant 
result reported in this article is the peculiar behavior 
of the quasiparticle density of states N{E) at the sur- 
face of the wire. Close to the surface, and for high total 
currents, the local density of states (LDOS) strongly de- 
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FIG. 3. Local density of states in the critical current con- 
figuration (/c) at different distances from the central axis for 
(a) a clean and (b) a dirty wire. The insets show the LDOS 
at the surface of the wire in (a) the clean wire, for currents 
I /lo = 0, 2.90, and 6.22 (critical value), and in (b) the dirty 
wire, for currents I/Io ~ 0,2.36, and 3.92 (critical value). 
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FIG. 4. Critical current /c as a function of radius R for 
three different degrees of disorder. Discrete points are the re- 
sult of the exact calculation. Lines are obtained from ansatz 
(^) in the main text for several values of F. The dash-dotted 
line gives the prediction from Silsbee's rule. Inset: Magnifi- 
cation of the small R quadratic behavior. 
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viates from its conventional BCS shape, to the point of 
losing the coherence peak and becoming gapless. This is 
shown in Fig. 3, where we plot the LDOS at different 
positions in the wire for the critical current configura- 
tion {I — Ic), as well as for different currents right at the 
surface (see insets). This strong distortion of the quasi- 
particle spectrum is caused by the large values which Vs 
acquires near the surface when the total current is high 
[^,0. On the contrary, near the center of the wire, 
where Vs is always small [see Fig. 2(b)], N{E) approaches 
the standard BCS form. In the dirty case, the fading 
of the gap at the surface is less pronounced, with the 
maximum value of iV(0) considerably smaller than in the 
clean case. This behavior occurs because disorder com- 
petes with Vs by tending to restore the gap in the LDOS 
1^,^. As a result, N{E) is less sensitive to the Vs, 
hence showing a weaker dependence on p and /. Unlike 
in the quasi-one-dimensional case jl^, a stable form of 
transport-induced gapless superconductivity is induced 
near the surface because of global stability. This strong 
transport dependence of N{E) at the boundary could be 
measured by performing a tunneling experiment on the 
surface of current-carrying wire, in the spirit of experi- 
ments made on superconducting films . We note that 
this is an intrinsic surface effect that will survive for wires 
of arbitrarily large radius. 

Fig. 4 shows the critical current Ic as a function of 
the radius R of the wire. A crossover from quadratic to 
linear behavior can be clearly appreciated as the radius 
increases beyond the penetration length. This crossover 
reflects the transition from quasi-one-dimensional to full 
three-dimensional behavior. Without the support of a 
fully self-consistent calculation like that presented here, 
an educated guess might be obtained from the assump- 
tions that j depends linear and locally on and that the 
critical current Ic{R) is reached when j{p ~ R) ^ ic(r), 
where jc(r) is the critical current density in a quasi-one- 
dimensional wire with disorder F |lo| , [l^ . Then we would 
obtain 

IciR)=27r\oMT)R^^^ (5) 
i-o[R/M) 

where Iq and Ii are modified Bessel functions of the first 
kind. Fig. 4 shows that ansatz (H) works well for small R 
(as expected, since then it is exact), but systematically 
underestimates Ic for i? > Ao in both clean and dirty 
wires. Silsbee's rule also yields values for Ic{R) which 
are too small | [T5[ . The extra current which the super- 
conductor accommodates as compared with the predic- 
tions of more simple treatments, can be traced to the 
non-monotonous radial dependence of the current den- 
sity, which allows for a bump near the surface. This 
result emphasizes the need for a fully self-consistent cal- 
culation in the formulation of a quantitative theory of 
the critical current. 

In summary, we have calculated the structure of the 



intrinsic Meissner effect in realistic wires. Our descrip- 
tion is based on a numerical procedure which entirely 
incorporates the non-local and non-linear aspects of the 
problem and which is exact within the context of the qua- 
siclassical formalism. The most relevant features are the 
existence of a peak in the radial profile of the current den- 
sity near the surface and the generation of a transport- 
induced gapless density of states at the surface, which 
could be measured in a tunneling experiment. These 
properties are displayed by stable high-current configu- 
rations in both clean and dirty wires. For wires of large 
radius, our calculation predicts values for the critical cur- 
rent greater than those obtained from phenonienological 
models. 
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